Universal and non-universal properties of wave chaotic scattering systems 
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The application of random matrix theory to scattering requires introduction of system-specific 
information. This paper shows that the average impedance matrix, which characterizes such system- 
specific properties, can be semiclassically calculated in terms of ray trajectories between ports. The- 
oretical predictions are compared with experimental results for a microwave billiard, demonstrating 
that the theory successfully uncovered universal statistics of wave-chaotic scattering systems. 



Wave systems appear in diverse branches of physics, 
such as quantum mechanics, electromagnetics and acous- 
tics. However, solving the wave equations can be difficult, 
particularly in the short wavelength limit Further- 
more, even if exact solutions were feasible, there may be 
uncertainties in the locations of boundaries or in param- 
eters specifying the system, and the desired wave quanti- 
ties can be extremely sensitive to such uncertainties when 
the wavelength is short. Thus, rather than seeking solu- 
tions for specific systems, it is often convenient to create 
statistical models which reproduce generic properties of 
the system This is the motivation for the applica- 
tion of random matrix theory to wave-chaotic systems, 
in which it is conjectured that useful statistics can be ob- 
tained by replacing the exact Hamiltonian or scattering 
matrices by random matrices drawn from an appropriate 
ensemble. Here, by wave-chaotic we mean that, in the 
small wavelength limit, the behavior of the wave system 
is described by ray orbit trajectories that are chaotic 3 . 
Although such formulations cannot predict details of any 
particular wave system, they do predict the distribution 
of properties in an ensemble of related wave-chaotic sys- 
tems. Random matrix theory is also hypothesized to pre- 
dict the statistical properties of a single wave-chaotic sys- 
tem evaluated at different frequencies (in, e.g., the cases 
of acoustics or electromagnetics) or energies (in the case 
of quantum mechanics). See Refs. 0, [y| for reviews of 
the theory, history, and the wide range of applications of 
random matrix theory. 

In this paper, random matrix theory is applied to 
model the scattering behavior of an ensemble of wave- 
chaotic systems coupled to the outside world through 
a single scattering channel (the generalization to larger 
numbers of scattering channels is straightforward [7|). 
Such scattering systems have been studied extensively, 
with most work focusing on the scattering parameter 
S d, [13, El, 113 1 which is the ratio between the re- 
flected waves and the incident waves in the channel. Here 
we consider ensembles of systems whose distribution of 
scattering parameters are well-described by the so-called 
Poisson kernel [1, [^, [13, [l^ . The Poisson kernel charac- 
terizes the probability density for observing a particular 
scattering parameter S in terms of the average scattering 
parameter S. It represents contributions to the scatter- 
ing behavior from elements of the system which are not 



random, such as the prompt reflection from the inter- 
face between the scattering channel and the chaotic sys- 
tem. In addition, rays within the scattering region which 
return to the scattering channel without ergodically ex- 
ploring the chaotic dynamics also affect S [ill . [T^ . [Tsj . 
The ability to determine S from first principles, thus in- 
corporating all non- universal effects, would dramatically 
improve our understanding and ability to uncover univer- 
sal fluctuations in a host of wave phenomena. Because S 
is the only parameter in the Poisson kernel, methods for 
finding it are of interest. Even though S can be estimated 
from experimental ensemble data, predicting it from first 
principles has so far not been addressed in general (al- 
though it has been done for the specific case of quantum 
graphs [l6j). In this paper we show how to semiclassi- 
cally obtain S, and we experimentally verify the accuracy 
and utility of our technique. 

A quantity equivalent to the scattering parameter S 
is the impedance, Z = Zo{l -f S)/{1 — S), where Zq is 
the characteristic impedance of the scattering channel. 
Because non-universal contributions manifest themselves 
in Z as simple additive corrections, we use Z in much 
of our discussion [l3, [l^, [l^ ■ We note that impedance 
is a meaningful concept for all scattering wave systems. 
In linear electromagnetics, it is defined via Ohm's law as 
V = ZI, where V represents the phasor voltage difference 
across the attached transmission line (the system's port) 
and / denotes the phasor current flowing through the 
transmission line. In acoustics, the impedance is the ratio 
of the sound pressure to the fluid velocity. A quantum- 
mechanical quantit y c orresponding to impedance is the 
reaction matrix (sl. Il7||. In what follows, our discussion 
will use language appropriate to scattering from a mi- 
crowave cavity excited by a small antenna fed by a trans- 
mission line (the setting for our experiments). 

With the transformation to impedance, if S is dis- 
tributed according to the Poisson kernel, we find that 
the impedance can be represented as Q 



Z — i-^avg ^avgi^j 



(1) 



where in the lossless case Ravg and Xavg are the real and 
imaginary parts of the impedance Zavg based on the aver- 
age scattering parameter, where Zavg = Zo{l-\-S)/{l — S) 
and iS, (which we call the normalized impedance) is a 
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Lorentzian distributed random variable with median 
and width 1. With uniform loss (e.g., due to an imaginary 
part of a homogeneous dielectric constant in a microwave 
cavity), Ravg and Xavg are the analytic continuations of 
the real and imaginary parts of Zavg as k — > k + ik/{2Q), 
where Q » 1 is the quality factor of the closed system, 
and k denotes the wavenumber of a plane wave. The 
normalized impedance of the lossy system has a uni- 
versal distribution which is dependent only on the ratio 
k/(2QAk), where Ak is the mean spacing between modes 

We find that Z„, 
classical limit 
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-^avg can be evaluated directly in the semi- 
Til as 



^avg 



Zr + Rr 




(2) 



where Rr, the radiation resistance, and Xr, the radia- 
tion reactance, are the real and imaginary parts of the 
radiation impedance Zr , which represents the impedance 
the system would have if all the energy which coupled 
into the system was absorbed before coupling back out 
, [l^, (...) indicates a suitable ensemble averaging (to 
be discussed further), b{l) is an index over all classical 
trajectories which leave the port and return to the port 
location with total path length I, Ci,(j_) is a function of 
the stability of the trajectory indexed by h{l), and Af,(;) 
is the action along the trajectory b{l) f3|. includes 
the initial phase shift and a geometrical factor that takes 
account of the spreading of the ray tube along its path. 
Here it has been assumed that the port radiates isotrop- 
ically from a location far from the two-dimensional cav- 
ity boundaries. These parameters can all be determined 
from the geometry of the scatterer and location of the 
port 0. 

The purpose of this letter is to test the accuracy and 
usefulness of Eq. 0. In practice, we take account of a fi- 
nite number of ray trajectories according to their length 
I. Therefore, in Eq. ^ we replace the summation by 
which signifies that the sum is now over all trajec- 

b(l),L 

tories b{l) with lengths up to L, I < L. 

In order to test our theoretical prediction, experimen- 
tal tests are carried out on a quasi-two-dimensional mi- 
crowave cavity with a single port [1, [l^j (see Fig. [U 
inset), where the length of Wall A is 21.6 cm, the length 
of Wall B is 43.2 cm, the distance of the port to the near- 
est wall (Wall D) is 7.5 cm, and the height of the cavity 
is 0.8 cm. For the frequency explored (6 to 18 GHz), 
higher-order vertical modes are cutoff so that the waves 
are quasi-two-dimensional. Furthermore, the correspond- 
ing wavelengths (1.7 cm to 5.0 cm) can be regarded as 
in the semiclassical regime, and the cavity shape yields 
chaotic ray trajectories. We excite the cavity by means 
of a single coaxial probe whose exposed inner conductor 
extends from the top plate and almost makes electrical 
contact with the bottom plate of the cavity 0]. 

The radiation impedance Zr in Eq. ^ is measured by 
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FIG. 1: (Color online) Plot of the impedance from the average 
of 100 cavity realizations (Z), versus frequency from 6 to 8 
GHz with perturbers inside the cavity. Shown are the real 
(three upper curves) and the imaginary parts (three lower 
curves) of the impedance for the theory {zi^l with L = 200 
cm, blue dashed) and the experiment {{Z), red solid), as well 
as the measured radiation impedance of the port {Zr, black 
thick). Inset: The wave chaotic two-dimensional cavity with 
perturbers and a single port. 



covering the four side walls of the cavity with microwave 
absorbers. Normalizing the measured impedance with 
the radiation impedance has been used to remove the 
non-universal properties due to the coupling of the port 
and the cavity [1, [2l|. Here we further consider the 
non-universal properties due to short ray trajectories by 
adding the summation term in Eq. ([2]). 

To verify that Eq. ([2]) describes non-universal char- 
acteristics of wave-chaotic systems, we first proceed to 
determine universal statistics by applying the ensemble 
average. Two irregular-shaped pieces of metal (with the 
maximum diameters 7.9 cm and 9.5 cm) are added as per- 
turbers in the wave-chaotic system that is shown in the 
inset of Fig. [1] where the circular dot shows the port and 
the two starlike objects represent the perturbers. The 
locations of the two perturbers inside the cavity are sys- 
tematically changed to produce a set of 100 realizations 
for the ensemble [20|. Typically, the shifts of resonances 
between two realizations are about one mean level spac- 
ing. The scattering parameter S is measured from 6 
to 18 GHz, covering roughly 1070 modes of the cavity. 
After the ensemble average, longer ray trajectories have 
higher probability of being blocked by the two perturbers, 
and therefore, the main non-universal contributions are 
due to short ray trajectories. We compare the measured 
ensemble averaged impedance (Z) and the theoretical 

impedance zi^l that is calculated from Eq. ^ with the 
sum up to the maximum trajectory length L — 200 cm, 
corresponding to a total of 584 trajectories. We multiply 
each term in the sum, Eq. (HI), by a weight equal to the 
fraction of perturbation configurations that are not inter- 
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FIG. 2: (Color online) The average x'^ of distributions of 
the phase of the scattering parameter ips on a semi-log 
scale, where the scattering parameters are calculated from 
impedance normalized with only the radiation impedance 
(green) and with ray trajectories according to the maximum 
trajectory length from L = 50 cm (red) up to L = 200 cm 
(blue), versus frequency window sizes from 0.1 to 4.0 GHz. 

cepted by the trajectory corresponding to that term. The 
result is shown in Fig. [T] where the three upper curves 
are the real part of the impedance (resistance) and the 
three lower curves are the imaginary part (reactance). 
The experiment curves (red solid) follow the trend of the 
radiation impedance curves (black thick) , and the theory 
curves (blue dashed) reproduce most of the fluctuations 
in the experiment curves with only a modest number of 
trajectories. The good agreement between the measured 
data and the theoretical prediction verifies that the new 
theory, Eq. ([J), predicts the non- universal features em- 
bodied in the ensemble averaged impedance. We believe 
that the differences between the measured data and the 
theory are due to the finite number of realizations in the 
ensemble and diffraction effects that are not taken into 
account in the theory. 

We further test our theory by consideration of the 
statistics of the scattering parameter for an ensemble of 
perturbation configurations and an ensemble of frequen- 
cies. Random matrix theory predicts that, after all non- 
universal effects have been removed, the resulting scat- 
tering parameter, which we denote as s = {i^ — l)/(iC + 
1) = |s|e''^% should be distributed uniformly in phase 
in < (fs < 27r, and this result is independent of loss, 
frequency, and mean level spacing 0, [l^ . The previous 
work of Refs. [1, [l^, [2l| removed the non-universal prop- 
erties by performing normalization with only the radia- 
tion impedance, as = {Z — iXpi)/ Rf>. Here we add ray 
trajectory corrections based on the maximum trajectory 
length L, 

z^(^) = (Z-zX(^))/i?(t), (3) 
and use the test to evaluate how uniform the result- 



FIG. 3: (Color online) Plot of the smoothed impedance versus 
frequency from 6 to 8 GHz. Shown are the real (three upper 
curves) and the imaginary part (three lower curves) of the 
smoothed impedance for the theory {Z^'"^ with L = 200 cm, 
blue dashed) and the experiment (red solid), as well as the 
measured radiation impedance of the port {Zr, black thick). 



ing phase distributions are. Ravg and Xavg are the an- 
alytic continuations of the real and imaginary parts of 
Za^g as A; — + fc -|- ik/(2Q) in the experimental case with 
loss. Experimental distributions of the phase ips of s are 
calculated from 100 realizations and different frequency 

1 ^ 

windows, x" ~ i — r / ("-i ~ {"'^ilf' nieasures the devia- 

tion between the experimental distributions of '-ps and a 
perfectly uniform distribution, where rij is the number of 
elements in the i*'' bin in the histogram (with ten bins, 
iV = 10) of the probability of the phase of the scattering 
parameter P{(ps), and (rii) is the average of rii over i. A 
smaller value means the experimental data are closer 
to the theoretical prediction. 

Fig. [2] shows the averaged evaluated over the spec- 
tral range from 6 to 18 GHz. The results indicate that 
the distributions of the measured data are systematically 
more uniform as more ray trajectories are taken into ac- 
count in the impedance normalization (Eq. (O). The 
improvement is dramatic after including just a few short 
ray trajectories {L 50 cm, 7 trajectories) and satu- 
rates beyond L = 100 cm (36 trajectories). The periodic 
wiggles represent the effects of the strongest remaining 
trajectory not taken into account in the theory. Thus, 
we see that non-universal effects of ray trajectories in 
the ensemble of wave-chaotic systems can be efficiently 
removed by considering a few short ray trajectories or by 
increasing the window size for the frequency ensemble. 

In addition to experiments with ensemble averaging 
over perturber positions, we now examine the theory in 
the stringent case of just a single realization without scat- 
terers, and we use only a frequency ensemble. We con- 
sider frequency smoothed experimental data and com- 
pare it with the smoothed theoretical prediction. Fig. 
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|3] shows that the smoothed measured impedance Z (red 
sohd) agrees with the smoothed theoretical impedance 

Z(^) (blue dashed). Notice that Z^^] -> Z(^) because 
there is only a single realization now. The smoothing 
function is a Gaussian with standard deviation tr = 240 
MHz, which inserts an effective low-pass Gaussian filter 
on the trajectory length, thus limiting the terms in Eq. 
^ to those with a path length / < c/cr = f 25 cm. The 
results in Fig. [3] shows that the theory correctly captures 
systematic contributions from short trajectories. 

In conclusion, the non-universal effects of coupling and 
short ray trajectories on wave scattering in chaotic sys- 
tems are predicted by a newly developed theory [t*! and 
verified experimentally. This is accomplished through 



statistical tests of the scattering parameter and compar- 
isons of impedance in an ensemble of perturbed systems, 
as well as a single-realization wave-chaotic system. In 
particular, non-universal effects have been better consid- 
ered and removed from measured data to reveal underly- 
ing universally fluctuations in the scattering parameter. 
These results should be useful in many fields, such as 
nuclear scattering, atomic physics, quantum transport in 
condensed matter systems, electromagnetics, acoustics, 
geophysics, etc. 

We acknowledge seminal discussions with R. E. Prange 
and assistance from Michael Johnson. This work was 
supported by the Air Force Office of Scientific Research 
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